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1. Introduction 

In one-particle-inclusive deep inelastic scattering (DIS) a lepton scatters off a nucleon and 
one of the hadrons produced in the collision is detected. In the one-photon exchange ap- 
proximation, the lepton-nucleon interaction proceeds via a photon of virtuality Q. The 
cross section depends in particular on the azimuthal angle of the final state hadron about 
the virtual photon well as on the azimuthal angle of the target polarization. In 

the kinematic region where the transverse momentum of the outgoing hadron is low com- 
pared to Q, the cross section can be expressed in terms of transverse-momentum-dependent 
parton distribution functions (PDFs) and fragmentation functions (FFs). These partonic 
functions are generalizations of the distribution and fragmentation functions appearing in 
standard collinear factorization. They are often referred to as unintegrated functions, as 
they are not integrated over the transverse momentum. 

The most complete treatment to date of one-particle-inclusive deep inelastic scattering 
at small transverse momentum remains the work of Mulders and Tangerman [jjl, com- 
plemented by Refs. ||2|, ^. In the last ten years, however, the subject has seen important 
theoretical and experimental progress. Initiated by the calculation of a nonvanishing Sivers 
effect in Ref . [Q] , unexpected developments arose with the correct treatment of Wilson lines 
in the definition of transverse-momentum-dependent PDFs and FFs |^, ^. In particular, 
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the fundamental tenet of universality of PDFs and FFs was revised ||7|, P]. New factor- 
ization proofs for the process under consideration here were put forward [|^, 11 1, updating 
past work ||l^. Some relations proposed in Ref. Q turned out to be invalid |l^, 14], and 



three new PDFs were discovered []15|, 16|. In the meanwhile, several experimental measure- 



ments of azimuthal asymmetries in semi-inclusive DIS were performed |17, 18, [T^, 2C, 21 



22, 23, 24, 25 



We consider it timely to present in a single, self-contained paper the results for one- 
particle-inclusive deep inelastic scattering at small transverse momentum, in particular 
including in the cross section all functions recently introduced. In Section ^ we recall the 
general form of the cross section for polarized semi-inclusive DIS and parameterize it in 
terms of suitable structure functions. In Section ^ we give the full parameterization of 
quark-quark and quark-gluon-quark correlation functions up to twist three and review the 
relations between these functions which are due to the QCD equations of motion. The 
structure functions for semi-inclusive DIS at small transverse momentum and twist-three 
accuracy are given in Section |^, and Section ^ contains our conclusions. The relation of 
the structure functions in the present paper with the parameterization in Ref. is given 
in Appendix and results for one-jet production in DIS are listed in Appendix ^. 

2. The cross section in terms of structure functions 

We consider the process 

m + N{P)^i{l') + h{Ph) + X, (2.1) 

where i denotes the beam lepton, the nucleon target, and h the produced hadron, and 
where four-momenta are given in parentheses. Throughout this paper we work in the one- 
photon exchange approximation and neglect the lepton mass. We denote by M and Mh 
the respective masses of the nucleon and of the hadron h. As usual we define q = l — l' and 
= —q^ and introduce the variables 

P-q P-Ph 2Mx 

Throughout this section we work in the target rest frame. Following the Trento conven- 
tions |8| we define the azimuthal angle (/>/i of the outgoing hadron by 



cos (ph = ; , sm (ph = , , (2.3) 

12 p2 //2 p2 

where lj_ = g^l^ and Pj^^ = Py^^^ are the transverse components of / and Ph with respect 
to the photon momentum. The tensors 

_ q^P'^ + P'^q" . 7^ fq'^q'' P^P' 



p.g(l+^2) + i + y Q2 ^2 )' i^-^i 



e^"^ = e^'^'P" ^^ "^^ (2.5) 

P-q^/l + j^ 
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Figure 1: Definition of azimuthal angles for semi-inclusive deep inelastic scattering in the target 
rest frame Ph± and S± are the transverse parts of Ph and S with respect to the photon 

momentum. 



have nonzero components = g'^ = —1 and = — Cj^ = 1 in the coordinate system of 
Fig. |l[ our convention for the totally antisymmetric tensor being e'^^^s _ ^ decompose 
the covariant spin vector S of the target as 



S'' = Si 



- qf" M'^ / {P ■ q) 



S-q M 



r 



S'^ = g^^S, (2.6) 



and define its azimuthal angle (j)s in analogy to (j)h in Eq. (|2.3[) , with Ph replaced by 5. 
Notice that the sign convention for the longitudinal spin component is such that the target 
spin is parallel to the virtual photon momentum for 5|| = —1. The helicity of the lepton 
beam is denoted by Ag. We consider the case where the detected hadron h has spin zero 
or where its polarization is not measured. 

Assuming single photon exchange, the lepton-hadron cross section can be expressed in 



a model-independent way by a set of structure functions, see e.g. Refs. |2^, We use 

here a modified version of the notation in Ref. |27|, see App. |A|, and write^ 



d(j 



dx dy dip dz dcph dP^j_ 




xyQ^ 2(1 



Fuu,T + eFuu,L + \/2e(l + e; 



cos (phF^"^^'' 



ecos{2(Ph) F;j°^"t"^ + Xe \/2e(l-6) 



sin (phF^j^^'^'' 



V^eil + e) sin F'Jl'^'^ + e sm{2M F'Jff"^ 



^The polarizations Sl and St in have been renamed to and \S±\ here. This is to avoid a clash 
of notation with Section ^ where subscripts L and T refer to a different z-axis than in Fig. 
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5'iiAr 



+ \Si 



Vl-e^ Fll + V2e(l-e) cos <Ph F'j^l '^^ 



,sm(3(f>h-4's) 



+ 15, |A, 



a/T^ cos(</);, - <^s) Fp^'^-^'^ + V2e(l-6) cos ./.^ 



icos 05 
LT 



LT 



(2.7) 



where a is the fine structure constant and the structure functions on the r.h.s. depend 
on X, Q^, z and F^_|_- The angle ijj is the azimuthal angle of I' around the lepton beam 
axis with respect to an arbitrary fixed direction, which in case of a transversely polarized 
target we choose to be the direction of S. The corresponding relation between and (j)s 
is given in Ref. |27|; in deep inelastic kinematics one has dil) ~ dcfis- The first and second 
subscript of the above structure functions indicate the respective polarization of beam and 
target, whereas the third subscript in F^^^, Fjju l and F^^^p^ '^^\ F^^'^ '^^^ specifies 
the polarization of the virtual photon. Note that longitudinal or transverse target polar- 
ization refer to the photon direction here. The conversion to the experimentally relevant 
longitudinal or transverse polarization w.r.t. the lepton beam direction is straightforward 
and given in [^]. The ratio e of longitudinal and transverse photon flux in i \i.7\) is given 

by 



e 



(2.8) 



so that the depolarization factors can be written as 

,.2 1 



y 



2(l-e) 1 + 72 
-2 1 



y 



2 l-e 



■ £ 



2(1-6) 



2(1-6) 



^2e{l + e) 
\/2e(l-e) 



2(1-6) 



1 + 7' 
1 

1+^ 
1 

1 

7it7 



{i-y-hV) 



{^-y + W). 
(i-y), 



(2 -y)Jl-y- l^^y^ « (2 - y) TT^, 



y\l'^-y- iT^y^ 



y{^-\y) 



yV^-y^ 
y{^-h)- 



(2.9) 
(2.10) 
(2.11) 
(2.12) 
(2.13) 
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Integration of Eq. (^]^) over the transverse momentum Ph± of the outgoing hadron gives 
the semi-inclusive deep inelastic scattering cross section 



da _ 2a^ 

dy d^/' dz 2 i 

+ |5x| v'2e(l + e) sin ^/.^ F^'^^^ + |5x|Ae V2e(l-e) cos ./.^ F^^).^ | , (2.14) 
where the structure functions on the r.h.s. are integrated versions of the previous ones, i.e. 

Fuu,Ti^, z, Q^) = j (fPh± Fuu,t{x, z, pI^,Q^) (2.15) 

and similarly for the other functions. 

We can finally connect the semi-inclusive structure functions to those for inclusive deep 
inelastic scattering. With the energies = {P-Ph)/M and v = {P-q)/M of the detected 
hadron and the virtual photon, we have 

da{iN^£hX) ^sc^-f^r^ da{lN ^ ihX) u + M da{m ^ IX) 



j dz dx dy dip v \ ^ ^ 



^ ^ dEh dx dy dij) v dx dy d'\\) 

(2.16) 

where we have summed over all hadrons in the final state, whose total energy is -|- M. 
Using {v -|- M)/v = 1 -|- 7^/(2x) we have 

-\Si_\\e^/2I{i^ cos(t)s2x^{gi+g2)^, (2.17) 

where 

^jdzzFuu,T{x,z,Q^)=2xFi{x,Q^) = Ft{x,Q'^), (2.18) 

h 

I dzzFuuA^^z^Q^) = i'^+7^)F2{x,Q^) -2xF^{x,Q^) = Fl{x,Q^), (2.19) 

h 

I dzzFLLix,z,Q^)=2x{giix,Q^)-^^g2ix,Q^)), (2.20) 
h •' 

j dzzFl^^'{x,z,Q^) = -2x^ {gi{x,Q^) + g2ix,Q^)) (2.21) 
h •' 

in terms of the conventional deep inelastic structure functions. For the relation with the 
more common expression for target polarization along or transverse to the lepton beam 
direction see Refs. [31, 27]. Finally, time-reversal invariance requires (see, e.g., Ref. [^ ]) 

Y j dzzFi;^^^x,z,Q^) =0. (2.22) 
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3. Transverse-momentum dependent distribution and fragmentation func- 
tions 



3.1 Light-cone coordinates 

Manipulations with parton distribution and fragmentation functions are conveniently done 
using light-cone coordinates. For an arbitrary four-vector v we write = {v^ it v^)/\/^ 
and vt = {v^jv"^) in a specified reference frame and give all components as [v~ ,vt]- 
We will use the transverse tensors g^^ and e^'^, whose only nonzero components are = 
= — 1 and = —e'^ = 1. The light-cone decomposition of a vector can be written in 
a Lorentz covariant fashion using two light-like vectors n+ = [0,1, Ot] and n_ = [1,0, Or] 
and promoting vt to a four- vector vt = [0,0, vt]- One then has 

yt^ = y+n'^+v-n'^+v!^, (3.1) 

where = v-ri-, v~ = v- and vt ■ ?^+ = vt ■ ?^- = 0. We further have 

gf = g^f^ - nlnt - n'^.n%, ef = e"'^"" n+pn_,. (3.2) 

Note that scalar products with transverse four-vectors are in Minkowski space, so that 
Vt ■ wt = —Vt ■ wt- 

For the discussion of distribution functions we will choose light-cone coordinates such 
that P has no transverse component, i.e. 

The spin vector of the target can then be decomposed as 

5^ = 5. (^-"'^<;/^-"^^"- +5^ (3.4) 

which implies Sl = M (S ■ n^)/{P ■ n^). Similarly, for the discussion of fragmentation 
functions, we will assume a coordinate choice with 

Pi: = PHn'- + §=<- (3.5) 

3.2 Calculation of the hadronic tensor 

We consider semi-inclusive DIS in the kinematical limit where becomes large while x, z 
and P^j^ remain fixed, and will perform an expansion in powers of 1/Q. For the calculation 
we use a frame where both ( ^^ ) and ( ^^ ) are satisfied, and where xP~^ = P^ I z = Q/y/2. 
Notice that this differs from the choice in Section |2|, where the transverse direction was 
defined with respect to the momenta of the target and the virtual photon, instead of the 
momenta of the target and the produced hadron. Details on the relation between the two 



choices can be found in [|l], 33|. In particular, Sl and St defined by ( |3.4| ) with (3^) and 
( |3.5| ) differ from S\\ and S^ in ( |2.6D by terms of order l/Q"^ and 1/Q, respectively. 
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The leptoproduction cross section can be expressed as the contraction of a hadronic 
and a leptonic tensor, 

dxdydi,dzd<PhdPl^ ^ 8^ '^^^'^^ ' ^^-^^ 
where the leptonic tensor is given by 

V = 2 (/^C + l'^lu-l- I'a^u) + 2iAe e^,p„ m'". (3.7) 
The hadronic tensor is defined as 

2MW^-^ = j^Yl I ^S('\q + P-Px-Ph){PmO)\h,X){h,X\rmP), (3.8) 

where J^{S,) is the electromagnetic current divided by the elementary charge and a sum is 
implied over the polarizations of all hadrons in the final state. 

The calculations in this paper are based on the factorization of the cross section into 
a hard photon-quark scattering process and nonperturbative functions describing the dis- 
tribution of quarks in the target or the fragmentation of a quark into the observed hadron. 
We limit ourselves to the leading and first subleading term in the 1/Q expansion of the 
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cross section and to graphs with the hard scattering at tree level. Loops can then only oc- 
cur as shown in Fig. ||b, c, d, with gluons as external legs of the nonperturbative functions. 
The corresponding expression of the hadronic tensor is |l|, |3^ ] 

2MW^"' = 2z^el I (fpT cPkT 5'^{pt + Qt - ^r) Ti4^''{x,pTh^A''{z, krh" 

(3.9) 

with corrections of order l/Q^ , where the sum runs over the quark and antiquark flavors 
a, and Ca denotes the fractional charge of the struck quark or antiquark. In the next 
subsections we discuss in detail the correlation functions $ for quark distributions, A for 
quark fragmentation, and their analogs and A^^ with an additional gluon leg. The first, 
second and third term in Eq. (|3.9D respectively correspond to the graphs in Fig. |2|a, b and 
c, with gluons having transverse polarization. The analogs of Fig. ^ and c with the gluon 
on the other side of the final-state cut correspond to the "h.c." terms in Eq. ( |3.9| ). 

The Wilson lines needed in the color gauge invariant soft functions come from graphs 
with additional gluons exchanged between the hard scattering and either the distribution or 
the fragmentation function (as in Fig. ^b,c,d). The corresponding gluons have polarization 
vectors proportional to in the first and to n_ in the second case, except for contributions 
from the gluon potential at infinity. The importance of the latter has been discussed in 
|5|, P , and a detailed derivation of the Wilson lines appearing in semi- inclusive DIS was given 
in Ref. [^] for the leading terms in the 1/Q expansion. For the contributions subleading 
in 1/Q only the cross section integrated over Ph± has been analyzed in the same reference. 

Going beyond the tree graphs just discussed requires modifications in the factorization 
formula ( |3.9D . In particular, radiative corrections involving low-frequency gluons introduce 
Sudakov logarithms and so-called soft factors. A proof of factorization to all orders in ag 
for the (similar but simpler) case of two-hadron production in e~^e~ collisions was given 
long ago |ll2[| . Recent work on all-order factorization in semi-inclusive DIS can be found 



in Refs. 11] and in Ref. Whether and how the tree-level factorization used in the 



present paper extends to subleading level in 1/Q is presently not known. 
3.3 The quark-quark correlators 

The quark-quark distribution correlation function is defined as^ 



(3.10) 



with p'^ = xP~^, where here and in the following we omit the flavor index a. The corre- 
sponding correlator for antiquarks is obtained by replacing the quark fleld by its transform 

^To be precise, one should distinguish the momentum fraction x in the definition of distribution functions 



from the Bjorken variable defined in Eq. (2.2). They coincide however in the process we consider, so that 
we drop this distinction for simplicity. An analogous remark holds for the argument z in the fragmentation 
functions below. 
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under charge conjugation, see Ref. 
lines) 



Q. In the correlator ( 3.10| ) we have gauge links (Wilson 



U. 



{0,+oo) 



OO 



(+00,0 y^T,^T^^')^'^'i^' ^C^t)- 
h~;cj') indicates a Wilson line running along the minus direction from 



(3.11) 
(3.12) 



Here U^^ia 
[a~,0+,CT] to [6~,0+,ct], while U'^{aT,bT;c 

ay] to [c~, 0"^ 



transverse direction from [c , O'^ 

W~{a-,b-:CT) 



V exp 



U'^{aT,bT;c )=Pexp 



) indicates a Wilson line running in the 
for], i.e. 

h- 



(3.13) 



(3.14) 



The correlator in Eq. ( |3.10| ) is the one appearing in semi- inclusive DIS. In different processes 
the structure of the gauge link can change 0, 34, 35, For instance, in Drell-Yan lepton 
pair production all occurrences of oo~ in the gauge links should be replaced by — oo~. 
In particular, this reverses the sign of all T-odd distribution functions appearing in the 
correlator (see below). In partonic processes with colored states in both the initial and 
final state, the gauge link contains contributions running to oo~ as well as — oo~, and 
T-odd terms differ by more than a simple sign change. Alternatively, it is possible to work 
always with the correlator for semi-inclusive DIS when convoluting T-odd functions with 
so-called gluonic-pole cross sections instead of the normal partonic cross sections |35, p7| . 
We note that beyond tree-level the Wilson lines in Eqs. ( 3.11 ) and (l3.12| ) lead to logarithmic 
divergences in the correlator ( 3.1C| ). These are due to gluons with vanishing momentum 
component along n_|_ and need to be regularized. The twist-two part of the correlator can 
be regularized in ways consistent with factorization to leading power in l/Q [|l2|, 3^, IC, 
It is presently not known how to extend this to the twist-three sector, see Ref. for an 
investigation of this case. 

A complete parameterization of the quark-quark correlation function has been given 
in Ref. |16|. Here we limit ourselves to the twist-three level, where we have 



+ hiT 



f: 



IT 



M 



+ hi 



[i>T, i 



2M 



2M 



+ 



M 
2P+ 



M 



PTpS'_ 



M 



Ta 



gTi5 $T + gs 75^ 



g 75 



per 

M 

per 

4 7pPTa 

M 



^-]75 r± [$T, [i^ 

2 2M 



(3.15) 
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The distribution functions on the r.h.s. depend on x and except for the functions with 
subscript s, where we use the shorthand notation [l| 

Pt ■ St 



9is{x,Pt) = SLgiL{x,PT) 



M 



■9it{x,Pt) 



(3.16) 



and so forth for the other functions. The first eight distributions of Eq. ( |3.15| ) are referred 
to as twist two, and the next 16 distributions are referred to as twist three, where we use 
the notion of "dynamical twist" as explained in Ref. [^]. The remaining eight functions 
of twist four have been omitted here and can be found in Ref. |^]. The 10 functions f^, 
hi, e^, Crp, e^, f!p, f^, f^, g^, h are T-odd |2|, |l^, i.e. they change sign under "naive time 
reversal" , which is defined as usual time reversal, but without the interchange of initial and 
final states. The functions g-^ [15|, and ||TH| exist because the direction of the Wilson 



line provides a vector independent of P and S for a Lorentz invariant decomposition of the 
correlation function ^{p), which is defined as in ( ^.lOD but with ^'^ integrated over instead 
of being set to zero [Q. The notation used in (|3.15| ) is consistent with Refs. |l|, |^, |l^, 16|, 



except for three points of discrepancy: (i) the sign of g is opposite to that in Ref. |15|, 
where the function was originally introduced, and consistent with Ref. |jl6|; (ii) the sign 



of the function is opposite to that in Ref. [16| and consistent with the other articles; 



(in) the function here is different from Ref. |16|, where it was first introduced, in order 
to maintain the symmetry with the other functions and to have simpler expressions in 
the following results. The relation between the two notations is the following (recall that 



Pt 



here 



Pt 
M2 



f±/ 
It 



Ref. [[| 



T 



here 



Ref. 16 



(3.17) 



The nomenclature of the distribution functions follows closely that of Ref. Q, sometimes 
referred to as "Amsterdam notation." We remark that a number of other notations exist 
for some of the distribution functions, see e.g. Refs. |41, 42, H^. In particular, transverse- 



momentum-dependent functions at leading twist have been widely discussed by Anselmino 



et al. 1 44, 45, 46 1. The connection between the notation in these papers and the one used 
here is discussed in App. C of Ref. 



We also list here the expressions for the traces of the correlator $(x,pt) from Ref. |16|. 



With = \ Tr[$r] we have 



$[7 ■''75] 



h 



€rp Ptp^t. 



Sl giL 



M 
Pt ■ St 

M 



IT ) 



■giT : 



$1' 



75J 



S^hi + SL^h 



IL 



a P 12 op 

PtPt ~ 2 Pt 9t 
M2 



Stp h 



IT 



cap 
^T PTp 

M 



n 5 



(3.18) 
(3.19) 

(3.20) 



M 



PTpS'_ 



M 



(3.21) 
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M 
P+ 

M 
P+ 



Pt-St 



ap 

ap Of e PTp j._ 



M 



Pt/'t 2Pt9t ^ naf±,PTfA 



M2 

brp 9T + SL— Ql 

a P 1 2 ctp 

PtPt ~ 2Pt9t q 



Tp9T 



M- 



ap 

PTp 



M 



M 

p+ 

M 



SrpJ)rj-i P'Y' Srp 

M 



Pt-St , 



(3.22) 



(3.23) 



(3.24) 

(3.25) 
(3.26) 



where a and (5 are restricted to be transverse indices. Here we made use of the combinations 



fT{x,PT) = M^^Pt) - 1^ fTix,PT), 

2 

9t{x,Pt) = 9t{x^Pt) - 1^9t{x,Pt), 



hi{x,p'^) = hiT{x,PT) 



Pt 
2M2 



^iT'(^'Pr)' 



(3.27) 
(3.28) 
(3.29) 



to separate off terms that vanish upon integration of the correlator over transverse mo- 
mentum due to rotational symmetry. The conversion between the expression in Eq. ( 3.23| ) 
and that in Eq. (19) of Ref. p| can be carried out using the identity 



Pt er^Tp = Pt ^tPtpStct + (pt ■ St) e^^'pTp , 



(3.30) 



which follows from the fact that there is no completely antisymmetric tensor of rank three 
in two dimensions. 

Integrating the correlator over the transverse momentum pT yields 



^>(x) = J (fpT Hx,pt) = ^ |/i 4+ + Sl9i 75 i+ + hi 



+ 1^ |e -iSLeL75 + h ^TSTpla + QT 75 9t 
where the functions on the r.h.s. depend only on x and are given by 



(3.31) 



h{x) = J d^PT fi{x,pl) 



(3.32) 
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and so forth for the other functions. We have retained here one common exception of 
notation, namely 

giix) = j (fpT giL{x,PT)- (3.33) 

Other notations are also in use for the leading-twist integrated functions, in particular 
f1 = q (unpolarized distribution function), gl = Aq (helicity distribution function), hi = 
6q = Axq (transversity distribution function). The T-odd functions vanish due to time- 
reversal invariance [llf 



J (fpT hix.pl) = 0, j (fpT eL{x,pl) =0, j cPpT h{x,pl) = 0. (3.34) 



We have kept them in Eq. ( 3.31 ) so that one can readily obtain the analogous fragmentation 



correlator, where such functions do not necessarily vanish. 
The fragmentation correlation function is defined as 



i-=o 

(3.35) 



with k = Py^ / z and the Wilson lines 



=^^^(oot,^t;+oo+)^/"+(W,^+;^t), (3.36) 



^("ot+oc) =^"^(0^'+°o^;0t)^/^(0t,oot;+oo+). (3.37) 

The notation U^+{a'^ ,b~^;CT) indicates a Wilson line running along the plus direction 
from [0~, ct] to [0^, 6+, ct], while U'^{aT, 6t; c+) indicates a gauge link running in the 
transverse direction from [0~, c^, ax] to [0~, c"^, bx]- The definition written above naturally 
applies for the correlation function appearing in e^e~ annihilation. For semi-inclusive DIS 
it seems more natural to replace all occurrences of +cxd^ in the gauge links by — oo^ p3|| . 
However, in Ref. Q it was shown that factorization can be derived in such a way that 
the fragmentation correlators in both semi- inclusive DIS and e"^e~ annihilation have gauge 
links pointing to +oo+. 

The fragmentation correlation function (for a spinless or an unpolarized hadron) can 
be parameterized as 

A(z,M = U^i^-+^^i"^^"' "^-^ 



2 y ' ' 2Mh 

where the functions on the r.h.s. depend on z and k?j.. To be complete, they should all 
carry also a flavor index, and the final hadron type should be specified. The correlation 
function A can be directly obtained from the correlation function <I> by changing^ 

n+^n_, er^-eT, P"^ ^ , M ^ Mh, x^l/z, (3.39) 



^The change of sign of the tensor er is due to the exchange n+ <-> n_ in its definition (p.2|). 



- 12 - 



and replacing the distribution functions with the corresponding fragmentation functions 
(/ is replaced with D and all other letters are capitalized). 

The correlator integrated over transverse momentum reads 

A(.) = z^/d^fc.A(z,M = I^i^ + ^{i^ + ^^i^^}, (3.40) 

where the functions on the r.h.s. are defined as 

Di{z) = j (fkTDi{z,kl) (3.41) 

and so forth for the other functions. The prefactor appears because Di{z,kT) is a 
probability density w.r.t. the transverse momentum k'j. = —zkx of the final-state hadron 
relative to the fragmenting quark. The fragmentation correlator for polarized spin-half 
hadrons is parameterized in analogy to Eqs. ( p. 15 ) and ( |3.31| ). As already remarked, in 



this case the functions Dt, El and H (the analogs fx, ei and h) do not vanish because 
\h,X) is an interacting state that does not transform into itself under time-reversal. Of 
course, it should be taken as an outgoing state in the fragmentation correlator. 

3.4 The quark-gluon-quark correlators 



We now examine the quark-gluon-quark distribution correlation functions |l33| , 47 



(3.42) 



which contain the covariant derivative iD^{^) = id^ + 9-^^ ■ Using U^j^^ = 
id~^ [^(+oo V'(0] can write 

<^+{x,pT)=xP+'^ix,pT) (3.43) 

for the plus-component of the correlator. For the transverse components we define a further 
correlator |Q 

ni^,PT) = n{x,PT)-p^^x,PT), (3.44) 

which is manifestly gauge invariant. It reduces to a correlator defined as in Eq. (|1|) with 
the covariant derivative iW^ replaced by gA^ if one has W^J^^ = 1, which is the case in a 
light-cone gauge = with suitable boundary conditions at light-cone infinity Q. Notice 
that ( |3.42| ) does not have the most general form of a twist-three correlation function since 
it depends on the kinematics of one instead of two partons (i.e., on a single momentum 
fraction and a single transverse momentum). Correspondingly, the covariant derivative 
is taken at the same space-time point as one of the quark fields. The tree- level result 
(|3.9| ) for semi-inclusive DIS can be expressed in terms of only the quark-quark correlator 
( 3.10D , the quark-gluon-quark correlator ( ^.42 ), and their fragmentation counterparts. This 



is not too surprising since the kinematics of our process defines a single plus-momentum 
fraction x, a single minus-momentum fraction z, and a single transverse momentum PhA_- 
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No further momentum fraction can for instance be constructed from the ratio -P^_|_/Q^, 
which is replaced with zero in the kinematical Umit we consider. 
The correlation function ( ^.44 ) can be decomposed as 



xM 



Ptp 



{Jis + ies)'~iTl5 + {h + ie) + (JiT - Ht) 



PTp ^Ta 

M 



M 



{97 



ap 



75) 



(3.45) 



where the index a is restricted to be transverse here and in the following equations. The 
functions on the r.h.s. depend on x and p^, except for the functions with subscript s, 
which are defined as in Eq. (|l6|). The last term in the curly brackets is irrelevant for the 
construction of the hadronic tensor of semi-inclusive DIS and has not been parameterized 
explicitly. The only relevant traces of the quark-gluon-quark correlator are 



2Mx 



Tr [!> 



Sl {hL + ieL) - ^^^^ (hr + ieT), 



Ij- Tr [1.^,(5^ + ^4^5)71 = f (/^ - ^g^) - S^p {It + ^9T) 



2Mx 



M 

Sl {R + ^^^^ ^I^t9t ^^^^ 5. ^ ^^X) ^ 



M 



where again we have used the combinations 

fT{x,PT) = fT{x,PT) 

9t{x,pt) = g'T{x,PT) 



Pt 
2M2 



fT{x,PT), 



Pt ~i-f 2 \ 



2M2 



(3.46) 
(3.47) 

(3.48) 

(3.49) 
(3.50) 



The above traces have been given already in Ref. ||T^ for the terms without transverse po- 
larization, whereas the terms with transverse polarization were partly discussed in Ref. ||l| 
(the functions ej^ and introduced in Ref. [16| were missing). 

Relations between correlation functions of different twist are provided by the equation 
of motion for the quark field 



[mo - m] m = b^iD'iO + + It iDaiO " m] m = 0, 



(3.51) 



where m is the quark mass. To make their general structure transparent we decompose 
the correlators into terms of definite twist, 



where the twist is indicated in the subscripts and the arguments {x,pt) are suppressed 
for ease of writing. One has = = and "P+'&^c 



0, where 



v. 



2l~l~^ and = 27^7^ projectors on good and bad Hght-cone compo- 

nents, respectively [40|. Projecting Eq. ( 3.51| ) on its good components one obtains for the 
correlators 



xM 7 ^>3 + M-^Tp^A 3 + |^T^2 - rn^2 



xM 7 $4 + M7j.p$^ 4 + i)T^ci- m^'i 



M 
P+ 



0, (3.53) 



where the term with has disappeared and the terms with and have been 
replaced using Eqs. ( |3.43| ) and ( ^.44 ). Multiplying this relation with one of the matrices 
r+ = {7^, 7^75, icr"^75}, which satisfy T'^V^ = r~^{l — V-) = r+, and taking the trace 
gives 



TrT^ 



m 
M 



■<i>9 



0, 



(3.54) 



where the terms multiplied by M/P^ in Eq. ( p. 53 ) have disappeared because the trace of 
Dirac matrices cannot produce a term that transforms like P"*" under boosts in the light- 
cone direction. Inserting the parameterizations ( p. 15 ) and (3.45) into ( p.54| ), one finds the 
following relations between T-even functions: 



771 



xe = 3;e + — /i, 
xf^ = xf^ + fu 



xgrp = xgrp +91T + Kt^ 
^9T = xgT-^9iT + j^h,, 



xgt = xgi + giL + 

xhL = xhL + j^h^L + —giL, 

7 , Pr , \ 
xhx = xhr -hi + + — giT, 



m 



xhrp = xhrp + hi + 



Pt 
2M2 



(3.55) 
(3.56) 
(3.57) 
(3.58) 
(3.59) 
(3.60) 
(3.61) 
(3.62) 
(3.63) 



These relations can be found in Ref. ||l|, App. C. Neglecting quark-gluon-quark correlators 
(often referred to as the Wandzura-Wilczek approximation) is equivalent to setting all 
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functions with a tilde to zero. For T-odd functions we have the following relations: 





(3.64) 


XCt = xct, 


(3.65) 


xerp — xep -+- — Jip, 

Ivl 


(3.66) 


2 


(3.67) 


x/t = ^/t + /iT) 


(3.68) 


xfT = xJt + /iT' 


(3.69) 




(3.70) 


xg =xg +j^hi, 


(3.71) 


xh = xh + ^-tJt hi~. 


(3.72) 



Most of these relations can be found in Refs. p8| , |49| ], and Eq. (3.71) can be inferred from 
Eq. (13) in Ref. |15|. Eqs. ( 3.66 ) and ( p. 68 ) have not been given before as they require 
the new functions introduced in Ref. |l^]. We emphasize that the constraints due to the 
equations of motion remain valid in the presence of the appropriate Wilson lines in the 
correlation functions. All that is required for the gauge link Z//{o,g) between the quark fields 
is the relation U(^Q^^)iD~^{^) = i9"'"Z//(o,$) leading to Eq. (|3.43|) . In contrast, the so-called 
Lorentz invariance relations used in earlier work are invalidated by the presence of the gauge 



links 1 14]. We remark that if quark-gluon-quark correlators are neglected, the time-reversal 
constraints ( 3.34 ) require that f d?Pp Pp fyp{x, pp) = and / (PppPphi{x,pp) = 0. 
The quark-gluon-quark fragmentation correlator analogous to is defined as 



A^),^.(z,fcT) 



-E 



(27r)3 



(3.73) 



The transverse correlator 



(3.74) 



can be decomposed as 



Mh 



+ {H + iE) i^- + ...{g^P-ie^P^,) 



i- 



(3.75) 
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Using the equation of motion for the quark field, the following relations can be established 
between the functions appearing in the above correlator and the functions in the quark- 
quark correlator ( 3.3^ ): 



E 

z 


E m 
= — + Di, 


(3.76) 


z 


= +Di, 

z 


(3.77) 


z 




(3.78) 


H 

z 




(3.79) 



4. Results for structure functions 



Inserting the parameterizations of the different correlators in the expression (3.9) of the 
hadronic tensor and using the equation-of-motion constraints just discussed, one can calcu- 
late the leptoproduction cross section for semi-inclusive DIS and project out the different 
structure functions appearing in Eq. ( |2.7D . To have a compact notation for the results, we 
introduce the unit vector h = Ph^./\Phl\ and the notation 

C[wfD]=xJ2^l [ d'PTd^kT6^^\pT-kT-PhJz)w{pT,kT)nx,pl)D^{z,kl), 

(4.1) 

where w{prp,krp) is an arbitrary function and the summation runs over quarks and anti- 
quarks. The expressions for the structure functions appearing in Eq. ( ^.71) are 



Fuu,T = C[hDi], 



Fuu,L = 0, 



Mh „ h-pT f ^Mh 



~ Q ^ 



psin 24>h _ n 



Mh ^ G^\h-pT 



fi- 



M 



T7 — [xeH, + —- 

Mh \ ^ M 

h-k^ f , rr±^Mh G^\^f^f .±n 
[h-krp] (h-pT) - krp-pj, . , 

WWh ''i^^i 



Mh ,^E\ 
-M^' -)\ 



FLL = C[giLDi], 



(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 

(4.8) 
(4.9) 
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pCOS 0^ _ '^M 



hk 



T 



Mh 



h prp 

M 



XQj Di + — — n-, J — 



(4.10) 



Fl 



f; 



sin(cf>h-(j>s) 
UT,T 

■sin{<l>h-<l>s) 
UT,L 

UT 

,sin(3<^h-0s) 
UT 



TTisin 4>s 



-,sin(2</>h-(/.s) 
UT 



M ■'^^ ^ 



h ■ kj 



2M 

krp * Pp 

2MMh 
2M 



xhTH, + —g,T — 



Q 



c 



2{h-pTf 



Pt 



2M2 



xf^Di 



2{h- krp^ (/l • prp) — krp • 

2MMh 



^1 M''^^ z 



M z J 
xhrH, + —g,T — 



xh^Hi 



-,cos{<j>^-4,s) 
LT 



^ LT 



+ 



hpj^ 



M 



F 



,cos{2<j>h-(j)s) 
LT 



2M ^ 

* Pp 

2MMh 
2M ^ 



giTDi 



xgrDi + hi 



Q 



xerHi — 
2{h-pTf 



M 
Mh 



z 



I I Mh , G-^ 

+ [-^THt + ^h-T- 



Pt 



2M2 



Mh 



E 



x9tDi + ^^iT— 



2 (h ■ krp) (h ■ prp) — krp ' prp 

2MMh 



xexHi 



Mh 



(4.11) 
(4.12) 
(4.13) 
(4.14) 



(4.15) 



(4.16) 
(4.17) 



(4.18) 



(4.19) 



Notice that distribution and fragmentation functions do not appear in a symmetric fashion 
in these expressions: there are only twist-three fragmentation functions with a tilde and 
only twist-three distribution functions without tilde. This asymmetry is not surprising 
because in Eq. ( |2.7| ) the structure functions themselves are introduced in an asymmetric 
way, with azimuthal angles referring to the axis given by the four-momenta of the target 
nucleon and the photon, rather than of the target nucleon and the detected hadron. 
Equations 



to ( 4.19 ) are a main result of this paper. A few comments concerning 
the comparison with the existing literature are in order here. First of all, it has to be 
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stressed that in much of the past hterature a different definition of the azimuthal angles 



has been used, whereas in the present work we adhere to the Trento conventions |28]. To 
compare with those papers, the signs of (j)h and of (f>s have to be reversed. The terms 
with the distribution functions and e^^ have not been given before. All leading-twist 
structure functions here are consistent with those given in Eqs. (36) and (37) of Ref. Q 
when only photon exchange is taken into consideration. The structure functions F^jj'^'^ 



and F^^^*^ in our Eqs. (|4.6|) and (4.7) correspond to Eqs. (16) and (25) in Ref. 



. The 

other six twist-three structure functions were partially given in the original work of Mul- 
ders and Tangerman [|^, but excluding T-odd distribution functions, assuming Gaussian 
transverse-momentum distributions, and without the contributions from the fragmentation 
function G^. 

The structure function F^^ '^^ is associated with the so-called Cahn effect [|^, |5l| . If 
one neglects the quark-gluon-quark functions D-^ in Eq. (4^) and /"*" in Eq. ( 3.56| ) as well 
as the T-odd distribution functions h and hj^, our result becomes 



hprp 

M 



(4.20) 



This coincides with the cos(j)h term calculated to order 1/Q in the parton model with 
intrinsic transverse momentum included in distribution and fragmentation functions, see 
e.g. Eqs. (32) and (33) in Ref. H. 

Let us briefly mention some experimental results and phenomenological analyses for 
the structure functions given above. For simplicity we do not distinguish between mea- 
surements of the structure functions and of the associated spin or angular asymmetries, 
which correspond to the ratio of the appropriate structure functions and Fuu^t + ^Fuu,L- 



1. 



Measurements of the cross-section components containing the structure function 
Fiju'^'^ have been reported in Refs. 54, 17, 2C]. A description of the coscph 
modulation by the Cahn effect alone has been given in Ref. fs^. The same analysis 
can be applied to the structure function F^°^ '^^ , leading to the results of Ref. 15^ 



pCOS 20;, 

tion 



contains the functions hf- (Boer-Mulders function 
). It has been measured in Refs. |2^ . 



]) and (Collins func- 



The structure function -F^j/ has been recently measured by the CLAS collabora- 



tion |5| 



4. The structure function F^^'^*' has been measured by HERMES [23|. The pre- 



cise extraction of this observable requires care because in experiments the target 
is polarized along the direction of the lepton beam and not of the virtual pho- 



ton [57, 58, 15, 27 1 . This implies that the longitudinal target-spin asymmetries mea- 
19|, |20|l receive contributions not only from F^^ '^'^ , but at the same 

''^"^ (see also the phenomenological 
In Ref. l2l the HERMES collaboration 



sured in Refs. ||1 
order in 1/Q also from F^^^ 
studies of Refs 



30, 61 



33|,|64|,^). 



has separated the different contributions to the experimental sin (/>/i asymmetry with 
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longitudinal target polarization and shown that F^'^ is dominant in the kinematics 



of the measurement. 



5. FjT^t!^ '^^^ contains the Sivers function p6| and has been recently measured for a pro- 



UT,T 



ton target at HERMES Q and for a deuteron target at COMPASS Extractions 
of the Sivers function from the experimental data were performed in Refs. |]67|, 68, B9| 



(see Ref. |7C] for a comparison of the various extractions). 



6. The structure function K 



,sin{</>h+0s) 



UT 



contains the transversity distribution function |41 



and the Collins function. As the previous structure function, it has been mea- 



sured by HERMES ||1 on the proton and by COMPASS g] on the deuteron. Phe- 



nomenological studies have been presented in Ref. |68|, where information about the 
Collins function was extracted, and in Ref. ||7l|, where constraints on the transversity 
distribution function were obtained by using additional information from a Collins 
asymmetry measured in e'^e~ annihilation |]72|| . 

Integration of Eqs. ( [4.2[ ) to ( [4.19 ) over the transverse momentum Ph± of the outgoing 
hadron leads to the following expressions for the integrated structure functions in Eq. ( 2.14| ): 



(4.21) 



Fuu,L = 0, 



(4.22) 



FLL = xY,elg1{x)Df{z), 



psin <j}s 



pcos </>5 
^ LT 



-X 



(4.23) 
(4.24) 



-X 



Y^el—i^ xgUx)Dl{z) + ^ h1{x) J . (4.25) 



Finally, the structure functions for totally inclusive DIS can be obtained from Eqs. ( 2.1(^ ) 
and (|1|) to (|]2l|). This gives the standard results |l|l 



2 

Fl = 0, 



9i = \^ei 9i{x), 

a 

91+92 = 1 ^elg!^{x), 



(4.26) 
(4.27) 

(4.28) 

(4.29) 



where given the accuracy of our calculation we have replaced gi — 7^52 by 51 in (4.28), and 
where we have used 



Y,fdzzDl{z) = l, J] /'dz^'^(z)=0, ^/'dzi?'^( 



z) = 0. 



(4.30) 
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The first relation is the well-known momentum sum rule for fragmentation functions. The 
second relation was already pointed out in Ref. f?^ . The sum rule for H follows from 
Eq. (4.24) and the time-reversal constraint ( 2.22| ). 



In App. y, we give results for one-jet production at low transverse momentum in DIS. 
5. Conclusions 

We have analyzed one-particle inclusive deep inelastic scattering off a polarized nucleon for 
low transverse momentum of the detected hadron, starting from a general decomposition 
of the cross section in terms of 18 structure functions, given in Eq. (2.7) and expressed 



in a helicity basis in App. Using tree-level factorization as discussed in Sec. |3.2| , the 
structure functions can be calculated up to subleading order in l/Q (twist three) using 
transverse-momentum-dependent quark-quark and quark-gluon-quark correlators, to sub- 
leading order, we need the parameterizations of the correlators up to twist three, involving 
24 parton distributions. In particular our treatment includes those twist three transverse 
momentum dependent functions that appear after the proper inclusion of Wilson lines 
in the quark-quark correlators. We also give the relations between the quark-quark and 
quark-gluon-quark correlators that follow from the QCD equations of motion. 

Using the parameterization of the correlators we eventually expressed the structure 
functions appearing in the cross section in terms of transverse-momentum-dependent par- 
ton distribution and fragmentation functions, see Eqs. ( |4.2| ) to ( [4.19| ). Several of these 
results were already present in the literature, but never collected in a single paper. The 
complete results for transversely polarized targets, Eqs. ( [4.11| ) to ( ^.19| ), including all T- 
odd distribution functions, are presented here for the first time. It is straightforward to 
generalize the present results to include the production of a polarized hadron. 
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A. Structure functions in a helicity basis 

Up to a kinematical factor, the structure functions we have introduced in Eq. ( |2.7] ) are 
simple combinations of cross sections or interference terms for the subprocess 7*A^ — > hX 
with definite helicities of the nucleon and the virtual photon. Let us define helicity structure 
functions 

ix q2 z Ph ) = ^'^^""'^ d + ^Y' (A 1) 



in terms of the j*p cross sections and interference terms introduced in Ref. [^7[, where j 
and i (n and m) are the helicities of the nucleon (virtual photon) in the amplitude and its 
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complex conjugate. We then have 
1 



Fuu,T 



rpcos (ph 



1 

7i 



Fuu,L = F++, 



(A.2) 



7?«i'^?^h - L Tm f F++ + F — 



(A.3) 



F^'^*^ = -^Im - F-o"), F^l^'^- = -ImF++, (A.4) 



F\ 



iCOS 



^Re(F+o+-F-o-), (A.5) 



,sin(0h-0s) 



f; 



,sin(0h+0s) 



(7r 



sin 



UT 



-ImF+-, 

Im F+r, 

-^Im F^ 
\/2 



+0 ' 



F, 



,j,sm(</>h-<^s) 
UT,L 

,sin(3</>h-0s) 



f; 



,sin(2</<h-0s) 



UT 



-ImFoV, 

-^ImF 
V2 



+0 ' 



(A.6) 



f: 



^COs(0h-0s) 

LT 

,cos{2if>h-<ps) 
LT 



ReF 



++ ' 



-pcos (ps 
LT 



(A.7) 



Comparing with our results Q to (p^ ) we find a number of simple patterns. The leading 
structure functions in the 1/Q expansion are those where the photon is transverse in both 
the amplitude and its conjugate, and the subleading structure functions correspond to the 
interference between transverse and longitudinal photon polarizations. The two structure 
functions Fjj^j ^ and F^™''^'" '^^^ involve only longitudinally polarized photons. In the 
kinematics we consider, they are of order 1/Q^ and thus beyond the accuracy to which 
we have calculated. We finally remark that the number of transverse momentum factors 
appearing in the different structure functions (|4.2D to ( [4.19| ) can be related to the mismatch 
between the helicity differences (m — i) and {n — j) using angular momentum conservation. 



see Ref. [g7|. 



B. Semi-inclusive jet production 

In this appendix, we take into consideration the process 

+ N{P) ^ £(/') + jet(F,-) + X 



(B.l) 
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in the kinematical limit of large at fixed x and -P^_|_. In the context of our tree-level 
calculation, we identify the jet with the quark scattered from the virtual photon. We then 
have z = \ and the cross section formula is identical to Eq. (^]^), except that it is not 
differential in z. Correspondingly, the structure functions do not depend on this variable. 
The structure functions for the process ( |B.1| ) can be obtained from those of one-particle 
inclusive DIS in Eqs. ( ^ ) to ( |4.19D by replacing Dx{z,k\) with 5(1 — 2;)(^(^)(fc2-), setting 



all other fragmentation functions to zero and integrating over z. This gives 



a 



2\P, 



a 

FLL = x^elglL{x,Pf^), 

a 

r^cos (hh 2 ^ -^i-L ±a/ t->2 \ 
a ^ 

j^s\n{cl>y,-4,s) _ _ 2 I-Pj-LI f^ar p2 x 



a 



^UT^^ _ X ^ — — xf^{x, Pjj_] 

a 

Pi 



B.2 

B.3 

B.4 

B.5 
B.6 

B.7 

B.8 

B.9 



psM2,,-,s) =xY,el ^ xf^'^ix^ Pf^), (B.IO 

a 

Fp'^^-^s^=xY,el^-^g^Mx,Ph)^ (B.ll 



^LT^' = E^f-^^aUx,Ph), (B.12 



^cos(2,,-,s) = _^ ^ xg^^{x,P^^), (B.13 



whereas the remaining 6 structure functions are zero. The results for the terms with indices 
UU, LL, and LT correspond to those in Ref. |^], Eqs. (119) to (121). The results for the 
terms with indices LU and UL correspond to those in Ref. ||l^. Integration of the cross 
section over Ph± leads to the results for inclusive DIS in Eqs. (|4.26| ) to ([4.29|) . Most 
terms vanish due to the angular integration, and F^^"^^ in Eq. ( [B.9D vanishes due to the 
time-reversal condition (3.34). 
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